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A general  empirical  expression  for  the  motion  of  a projectile  of 
mass  m and  velocity  V obeying  a seoond  order  drag  law  has  been  pro- 
posed that  has  the  form1 


F = mjjf  = - A'V2  - B’V  - C'  . (1) 

The  coefficients  A',  B',  and  C'  are  quantities  that  can  be  either 
evaluated  experimentally  from  penetration-time  data  or  calculated 
using  basic  quantities.  The  terms  on  the  right  side  of  Eq.  1 can  be 
considered  to  be  the  forces  contributing  to  three  components  of  energy 
deposition.  These  three  components  or  modes  of  energy  deposition  are: 

(1)  the  acceleration  of  target  material  adjacent  to  the  projectile, 

(2)  the  effect  of  frictional  forces,  and  (3)  the  cohesive  strength  of 
the  target  material  respectively.  A comprehensive  bibliography  of 
earlier  work  using  versions  of  Eq.  1 where  one  or  two  of  the  coef- 
ficients were  set  to  zero  is  given  by  Goldsmith2  on  p.298. 

Dienes  and  Welch ^ define  the  coefficients  in  terms  of  basic 
physical  quantities.  Their  definitions  are: 

A'  = J PC/  (!A) 

B'  = 6irrn  (IB) 

C'  = itt2o  (1C) 

where  p is  the  target  density  and  nr2  is  the  frontal  area  of  the 
projectile.  The  quantity  is  the  drag  coefficient,  n is  the 
dynamic  viscosity,  and  o isuthe  normal  stress. 

The  rate  of  work  by  the  projectile  as  it  penetrates  the  target 
is  equal  to  its  rate  of  loss  of  kinetic  energy  and  is  given  by 

F • dS  = d(j  mV2)  (2A) 


1’W.A.  Allen , E.B.  May  fields  and  N.L.  Morrisons  "Dynamics  of  a Projec- 
tile Penetrating  Sands"  Journal  of  Applied  Physics,  28,  1957. 

2 

'Weiner  Goldsmith,  "Impact,"  Eduard  Arnold  (Publishers)  Ltd.,  London, 
England,  1960. 

3 

' J.K.  Dienes  and  J.E.  Welch , "Parameter  Study  of  Rigid  Body  Pene- 
tration," Systems,  Science  and  Software  Report  No.  SSS-R-7 5-2521, 
1974. 
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This  leads  in  the  limit  to 

F = mV^r  = - A'V2  - B'V  - C'  . (2B) 

Eqs.  1 and  2B  are  solved  for  ?n  impacting  velocity  Vq  at  time  zero 

to  obtain  analytical  relations  for  the  velocity  and  penetration  as 
functions  of  time.  The  solutions  are  discussed  in  the  next  section. 

A computer  code  NEWT  was  written,  using  the  BASIC  programming 
language,  to  calculate  the  penetration  and  velocity  values  for  a 
given  impacting  velocity  and  a given  set  of  coefficients  at  times 
subsequent  to  impact.  The  energy  deposited  (both  total  and  the 
components  contributed  by  each  force  component)  along  the  projectile 
trajectory  are  also  calculated. 

A listing  of  the  code  is  given  in  Appendix  B and  a description 
of  the  code  that  includes  a flow  chart  is  given  in  Section  IV.  A 
description  of  the  code  input,  to  be  supplied  by  the  user,  is  given 
in  Section  V.  Sample  problems  that  demonstrate  the  differing  effect 
of  each  coefficient  on  the  projectile  penetration  are  discussed  in 
Appendix  A. 

The  code  has  been  run  successfully  on  the  WANG  2200  calculator. 
Results  and  identifying  captions  are  printed  by  the  calculator- 
controlled  typewriter.  Approximately  5000  words  of  memory  are 
needed  to  execute  the  program. 


II.  CALCULATIONAL  PROCEDURE 
Equation  1 may  be  expressed  in  integral  form  as 


V T 

J M?  'WTc  * - / 4r  = -T  • (3) 

V o 

o 

where  Vq  and  V are  the  velocities  at  time  zero  and  T respectively. 

The  constant  m has  been  absorbed  into  the  coefficients  yielding  the 
unprimed  set.  The  integrals  can  be  solved  in  closed  form  where  the 
solution  to  the  leftmost  integral  may  assume  one  of  the  three  forms 
dependent  upon  the  value  of  the  discriminator  R1  where 

R1  = B2  - 4 AC  . (4) 
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T = I *n  0 

1 R 2AV  ♦ B + R 


which  expands  to 


where 


1 /l  2AV  + B ♦ r\ 
" R Ml \ G 2AV  + B - R J 


2AV  + B + R 
o 

2AV  + B - R 
o 


The  term  G will  always  be  positive.  The  velocity  is  given 
explicitly  by 

„ R + B + GeRT(R  - B) 


2A(GeRT  - 1) 

That  relation,  as  expected,  will  give  a value  of  Vq  to  the 
velocity  when  T is  zero. 

The  solution  to  Eq.  3 for  negative  values  of  R1  is  given  by 

V 

rp  2 tan-1  2AV  + B ° 

R R 


(12A) 


(12B) 


which  expands  to 


2 -1  »».  * B , -1  2AV  * b\ 

r\,m  — e tan  — r — / 


The  velocity  is  given  explicitly  by 


5a  [ 


. 2AV  + B 

tan  (tan'1  2 |L  ) 


_R]  * 


The  solution  to  Eq.  3 for  a zero  value  of  R1  is  given  by 

T.  |V 

2AV  + B 


which  expands  to 


T = 2 [ 2AV  ♦ B " 2AV  + B ] 


The  velocity  is  given  explicitly  by 

4AV  - BT (2AV  ♦ B) 

V 0 0 

2A[2  ♦ T(2AVq  + B)] 
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The  preceding  solutions  are  valid  only  for  positive,  non-zero  values 
of  the  coefficient  A.  The  integral  in  Eq.  3 was  solved  for  zero  values 
of  A and  positive,  non -zero  values  of  B yielding 


v = l [e-BT(Bv0.C)  -c]  . 

Similarly,  Eq.  8 was  solved  for  the  penetration  yielding 


v * F*n 


BV  + C 

BV  + C 
o 


(20) 

(21) 


None  of  the  preceding  solutions  are  valid  if  both  A and  B are 
zero.  For  this  case  the  solution  to  Eq.  3 was  found  to  be 


V = V CT  . 
o 


The  solution  for  the  penetration  was 


(22) 


(23) 


III.  CODE  DESCRIPTION 

The  code  NEWT  was  designed  so  that  it  would  check  the  input 
values  of  A,  B,  and  C,  evaluate  the  discriminator  if  necessary, 
choose  the  proper  set  of  equations  to  use,  and  calculate  the  pene- 
tration, velocity,  and  energy  deposition  rates  at  the  chosen  times. 
The  input  and  the  calculated  results  are  then  printed  in  tabular 
form  along  with  their  captions. 

A flow  chart  of  NEWT  is  displayed  in  Figure  1.  The  important 
variables  and  arrays  are  listed  and  briefly  described  in  Table  I. 


-***  ■****&» 


Start 


Read  Input  Data 


Determine  Calculation  Mode 


A=0,  B>0,  C>0  Choose^^  A=0,  BsQ#  c>0 

I \Calculation  Mode^> 1 


A>0,  B^O,  C=0  Calculate  Velocity 


<C^More  — 
| No 

Calculate  Energy 
Disposition  Rates 

i , 

Print  Results 


Figure  1.  NEWT  Flow  Chart 
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Table  I.  Variables  and  Arrays 


i 


Variable 


Description 


A3 

B3 

C3 

A 

B 

C 

R 

M 

V 

o 

Q7 


FI 


N1 

T 

T1  (50) 
V 

VI (50) 
S 

SI (50) 
R1 


The  coefficient  A in  Equation  1 
The  coefficient  B in  Equation  1 
The  coefficient  C in  Equation  1 
A = A3  M/0. 5nR2 
B = B3  M/6ttR 
C = C3  M/ttR2 

Effective  impacting  radius  of 
Projectile 

Mass  of  Projectile 

Initial  Velocity  of  Projectile 

Calculation  mode  flag 


Q7  = 

001, 

A = 

B = 

0,  C > 0 

Q7  = 

010, 

A = 

C = 

0,  B > 0 

Q7  = 

011, 

A = 

C = 

0,  C > 0 

Q7  = 

100, 

A > 

o. 

B = C = 0 

Q7  = 

101, 

A > 

o. 

B = 0,  C > 0 

Q7  = 

111, 

A > 

o, 

B > 0,  C > 0 

Calculation  discriminator 
FI  = 1,  B2  -4 AC  > 0 

FI  = 2,  B2  -4AC  < 0 

FI  = 3,  B2  -4 AC  = 0 

Number  of  times 
Times 

Times  array 
Velocity 
Velocity  array 
Penetration 
Penetration  array 
R1  = /±(B2  - 4 AC) 
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IV.  INPUT 


■ 

i 

! 


The  user-supplied  input  needed  to  execute  the  program  is  listed 
in  Table  II.  It  is  listed  in  the  order  by  which  it  is  requested  by 
NEWT.  C6S  units  should  be  used  to  conform  to  captions  but  they  could 
be  easily  changed  to  conform  to  another  system.  It  should  be  noted 
that  the  programming  variables  A,  B,  and  C are  not  the  coefficients 
A,  B,  and  C.  The  coefficients  A,  B,  and  C cor-espond  to  the  programming 
variables  A3,  B3,  and  C3. 


Table  II.  Code 

Instruction 

1.  Title  = 

2.  A = 

3.  B = 

4.  C = 

5.  Radius  = 

6.  Mass  = 

7.  Initial  Velocity  * 

8.  No  Times  = 

9.  1 for  Read-in  Times  = 

10.  Times  Interval  = 

^11.  Next  Time  = 


•Generated  Queries 

Description 

Problem  Title 
A = pCD  (Eq.  1A) 

B = n(Eq.  IB) 

C = cr(Eq.  1C) 

Frontal  Area  = rrr2 
Mass  = (gm) 

VQ  = (cm/ sec) 

No  of  times  for  which  kinetic 

quantities  are  to  be  calculated 
F2  = 0 Calculates  times 
Time  Interval  (calculated  times  only) 
Next  Time  (Read-in  Times  only) 


V.  CONCLUSIONS 

This  study,  and  the  building  of  the  Computer  Code,  were  conducted 
at  the  Target  Assessments  Branch  (TAB),  Vulnerability/Lethality 
Division (VLD) , Ballistic  Research  Laboratories  (BRL)  as  part  of  an 
on-going  effort  to  calculate  the  penetration  of  human  tissue  by  bullets 
or  fragments  and  to  derive  the  resulting  incapacitation  probabilities.  ' 

However,  the  code  is  not  restricted  to  a tissue  media,  but  can  be  used 
to  calculate  the  penetration  of  any  media  that  is  described  by  a second 
order  drag  law.  A magnetic  tape  copy  of  the  code  can  be  obtained  by  • ?v 

sending  a blank  computer  tape  cassette  to  the  TAB. 
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Earlier  work  by  Bruchey*4,  in  performing  least  square  determination 
of  drag  coefficients  for  handgun  bullets  penetrating  tissue  simulant, 
has  provided  a validation  of  the  model  for  non- zero  values  of  the  A- 
and  C- coefficients  and  a zero  value  of  the  B- coefficients.  It  would 
therefore  be  expected  that  a better  fit  of  experimental  data  to  the 
model,  as  well  as  a better  description  of  the  energy  deposition  modes 
in  the  medium,  could  be  obtained  by  optimizing  all  three  coefficients 
instead  of  only  the  A- coefficients  as  was  done  in  the  preceding  work. 

A computer  code  for  providing  this  least  squares  fit  of  experimentally- 
determined  penetration- time  data  to  the  model  by  optimizing  all  three 
coefficients  is  being  developed  by  the  TAB.5 

The  degree  of  human  incapacitation  (however  this  term  is  ultimately 
quantized  and  defined)  will  probably  be  the  sum  of  the  products  formed 
by  the  volume  of  a particular  tissue  destroyed  and  its  vulnerability 
index.  The  volume  of  tissue  destroyed  is  approximated  by  the  volume 
of  the  maximum  temporary  wound  cavity  (TWC).  A future  effort  will  be 
made  to  construct  a wound  cavity  model  that  used  the  energy  deposition 
rates  to  predict  the  size  of  the  TWC.  It  is  expected  that  the  forth- 
coming study  and  the  resulting  Computer  Code  will  provide  a rapid  and 
inexpensive  means  of  calculating  the  volume  of  that  cavity. 


4 

' William  J. Bruchey,  Jr.,  "Effective  Drag  Coefficients  for  Various  Small 
Arms  Projectiles  in  Dense  Media,"  BEL  Memorandum  Report  2481 , May  1975. 
(AD  #B004154L) 

’William  B.  Beverly,  " LSP : A Least  Square  Program  for  Calculating  The 
Resistive  Force  Coefficients  of  a Projectile  Subject  to  a Second  Order 
Drag  Lau, " to  be  published  as  a BRL  Report. 
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APPENDIX  A 
SAMPLE  PROBLEMS 


Sample  values  A,  B,  and  C were  chosen  so  that  the  three  components 
of  the  opposing  force  would  be  equal  at  the  impacting  velocity  of  30,000 
cm/sec.  This  choice  of  input  values  permits  an  easy  analysis  of  the 
kinematic  results  for  demonstrating  the  differing  effects  of  the  three 
force  components.  The  penetration  and  velocity  versus  time  are  listed 
in  Table  A-I,  the  partial  differential  (in  space)  energy  losses  are 
listed  in  Table  A-II,  and  the  partial  cumulative  energy  losses  are 
listed  in  Table  A-III. 

i 

i 

Three  more  calculations  were  conducted  in  which  A,  B,  and  C in 
turn  were  assigned  a non-zero  value  (that  used  in  the  first  problem) 
while  the  remaining  two  coefficients  were  set  to  zero.  The  penetration 
and  velocity  versus  time  values  calculated  for  these  problems  are 
listed  in  Tables  A-IV  - A-VI.  The  penetration  and  velocity  versus  time 
for  all  test  problems  are  displayed  in  Figure  A-l. 


* 


Table  A-I.  Kinematic  Results  (A>0,  B>0,  OO) 


TEST  PROBLEM  1 


MASS**  8. 

000E+O0  RADII’S- 

6.000E-01  INI? 

VEL-  3.000F+0/ 

A-  1. 

, 200F+00  B— 

1.800E+03  O 

5 . 399F+00 

A3=  8. 

482E-02  B> 

2.544E+03  C3*= 

7.633E+07 

KINEMATIC  RESISTS 

TIME (SEC) 

PEN  (CM) 

VEL  (Cl  I/SEC) 

DEL  E/DEL  S 

O.OOOOE+OO 

O.OOOOE+OO 

3.0000E+04 

2.2901E+08 

1.0000E-05 

2.8883E-01 

2.7793E+04 

2.1258E+02 

2.0000F-05 

5.5f)38E-01 

2.5741E+04 

1.9805E+08 

3.0000E-05 

8.0412E-01 

2.3827E+04 

1.8512E+08 

4 . OOOOE-05 

1.0533E+00 

2.2O34E+04 

1. 735°EJ-00 

5.000012-05 

1.2451E+00 

2.O351E+04 

1.6825E+0" 

C.O0O0E-05 

1.4406E+00 

1.3766E+04 

1.53OC,e+0" 

7.00OO”-O5 

1.6208E+00 

1.7269E+04 

1.4557E+0". 

3.000?  ”-05 

1.78r.3E+oo 

1.5351E+04 

1. 37°°E+0r! 

9.000012-05 

1.938012+00 

1.4506E+04 

1.3110E+O" 

1.0000E-O4 

2.07G7E+0O 

1.3227E+04 

1.2484E+O0 

1.1000E-04 

2.2028E+00 

1.2008E+04 

1.1912E+09 

1.2000E-04 

2.3170E+00 

1.0843E+04 

1.X3O0F+0'’ 

1 . 3000E-04 

'’.4193E+00 

0.7287E+03 

l.°r’12E+08 

1.400012-04 

2.5117E+00 

8.6597E+03 

1.0473E+08 

1.5000F— ''4 

2.5932E+00 

7.632°iE+03 

I.007OE+O0, 

1.6000E-04 

2.6645E+00 

6.G446E+03 

9.69°lE+07 

1.7000E-04 

2.7262E+00 

5.6920E+03 

°.  3 57^+07 

1.30Q0E-04 . 

2.7735E+0 o 

4.7723E+03 

°.O413E+07 

1.9000E-04 

2.3217E+00 

3. "029E+03 

8.74°7F+07 

2.000012-04 

3.8562-2+00 

3.021612+03 

8.4801E+07 

2.1000E-04 

2.8822E+00 

2.1862E+03 

8.2306E+07 

2.2000E-04 

2 . 9000E+00 

1. 3749E+03 

7.OP''6E+r'7 

2 . 3000E-04 

2.O090E+00 

5.3577E+02 

7. 7857FX37 
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Table  A-II.  Differential  Energy  Deposition  (A>0,  B>0,  C>0) 


DIFFERENTIAL  ENERGY  LOSS (ERGS /CM) 


TIME (SEC) 

A3*V*V 

B3*V 

r 3 

O.OOOOE+OQ 

7.6340E+07 

7.6340E+07 

7.6337E-HJ7 

1.0000E-05 

6.5524E+07 

7.0726E+07 

7.8337E+07 

2 . OOOOE-05 

5.6207E+07 

6.5505E+07 

7.6337E+07 

3. 0000E-05 

4.8157E+07 

6.0633E+07 

7.6337E+07 

4. 0000E-05 

4.1184E+07 

5.6071E+07 

7.6337E+07 

5. OOOOE-05 

3.5131E+07 

5.1788E+07 

7.8337F-M)? 

6.O0OOE-O5 

2.9871E+07 

4.7753E+07 

7.6337E+07 

7. OOOOE-05 

2.5296E+07 

4.3944E+07 

7.6337E+07 

8. OOOOE-05 

2.1314E+07 

4.0338E+07 

7.6337E+07 

9. OOOOE-05 

1.7851E+07 

3.6915E+07 

7.6337E+07 

1.0000E-04 

1.4841E+07 

3.3660E+07 

7.0337E+07 

1.1000E-04 

1.2231E+07 

3.0557E+07 

7.6337E+07 

1.2000E-04 

9.9737E+06 

2.7593E+07 

7.A337E+07 

1.3009E-04 

8.0283E+06 

2.4756E+07 

7.6337E+07 

1.4000E-04 

6.3610E+06 

2.2036E+07 

7.6337F.+07 

1.5000E-04 

4.9418F+06 

1.9423E+07 

7.G337E+07 

1.G000E-04 

3.7450E+06 

1.6908E+07 

7.G337E+07 

1.7000E-04 

2.7482E+06 

1.4484E+07 

7.6337E+07 

1 . 8000E-04 

1.9318E+OG 

1.2144E+07 

7.6337E+07 

1.9000E-04 

1.2789E+06 

9.8809E+06 

7.G337E-K)7 

2.0000E-O4 

7.7447E+05 

7.6891E+06 

7.6337F+07 

2.1000E-04 

4.0544E+05 

5.5G34E+06 

7.0337E+07 

3.2000F-04  ' 

1.6034E+05 

3.4987E+0G 

7.6337E+07 

2.3000E-04 

2.9105E+04 

1.4906E+06 

7.6337E-W 

* 

« 
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Table  A-III.  Cumulative  Energy  Deposition  (A>0,  B>0,  C>0) 


CUMULATIVE  EIIERGY  LOSS  (ERGS) 

TIME (SEC) 

INT (A3*V*V) DS 

Iirr(B3*V)DS 

IUT(C3)DS 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

1.0000E-05 

1.6390E+08 

1.6991E+08 

1.7639E+OR 

3.0000E-05 

2.9418E+08 

3.1.r>70E+08 

3.3973E+03 

3 . 0000E-P5 

3.9760E+08 

4.4070E+03 

4.9100E+07 

4.0000E-05 

4 . 7951E+08 

5.4770E+08 

6.3106E+00 

5 . 0000E-O5 

5.4418E+08 

6.3910E+08 

7.6043E+03 

6.0000E-05 

5.9502E+08 

7.1695E+08 

3.7983E+08 

7.0000E-05 

6.3476E+08 

7.8301E+08 

9.8932E+08 

8.0000E-05 

6.6562E+08 

8.3882E+08 

1.090OE+09 

9.0000E-05 

6.8940E+00 

8.8571E+08 

1.1835E+09 

1.0000E-04 

.7.0752E+08 

9.2484E+08 

1.2632E+09 

l.lOOOE-04 

7.2118E+08 

9 . 5724E+08 

1.3452E+09 

1.2000E-04 

7.3133E+08 

9.8381E+08 

1.4150E+09 

1.3000E-04 

7.3873E+03 

1.0053E+09 

1.4778E+09 

1.4000F-04 

7.4402E+08 

1.0225E+09 

1.5339E+09 

1.5000E-04 

7.4770E+08 

1.0360E+09 

1.5836E+09 

1.6000E-04 

7 . 5018E+08 

1.0464E+09 

1.6272E+09 

1.7000E-04 

7.5178E+08 

1.0541E+09 

1.6649E+09 

1.8000E-04 

7.5276E+08 

1.0597E+OO 

1.69G8E+09 

1.9000E-04 

7.5332E+08 

1.0635E+09 

1.7232E+09 

2.0000E-04 

7.5360E+08 

1.0659E+09 

1.7443E+09 

2.1000E-04 

7.5372E+08 

1.0673E+09 

1.7602E+09 

2.2000E-04 

7.5376E+08 

1.0679E+09 

1.7710E+00 

2 . 3000E-04 

7.5377E+03 

1.0681E+09 

1.7770E+0O 
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Table  A- IV.  Kinematic  Results  (A>0, 

B=0 , C=0) 

TEST  PROBLEM  2 

MASS-  8 

.OOOE+OO  RADIUS- 

G.OOOE-Ol  INIT 

VEL-  3.000E+04 

A-  1 

.20QE+00  B- 

0. OOOE+OO  C- 

0. OOOE+OO 

A.3-  8.482E-02  B3- 

0. OOOE+OO  C3- 

0. OOOE+OO 

KINEMATIC  RESULTS 

: 

TIME  (SEC) 

pen  (at) 

VEL(CM/SEC) 

DEL  E/DEL  S 

O.OOOOE+OO 

O.OOOOE+OO 

3.0000E+04 

7.6340E+07 

1.0000E-05 

2.9624E-01 

2.9255E+04 

7.25°"E+07 

2.0000E-05 

5.8523E-01 

2.8547E+04 

6.9125E+07 

3.0000E-05 

8.6730E-01 

2.7872E+04 

6.5895E+07 

4 . OOOOE-05 

1.1427E+00 

2.7228E+04 

6.2P.86E+07 

5. 0000E-05 

1.4119E+00 

2.6613E+04 

6.O079E+07 

6. OOOOE-05 

1.6751E+O0 

2.6026E+04 

5.7456E+07 

7. 0000E-05 

1.9325E+00 

2.5464E+04 

5.5001E+07 

8. 0000E-05 

2.1845E+00 

2.4925E+04 

5.2699E+07 

9. 0000E-05 

2.43I1E+00 

2.4409E+04 

5.0540E+07 

1.0000E-04 

2.6727E+00 

2. 3914E+04 

4.8510E+07 

1.1000E-04 

2.9095E+00 

2.3439E+04 

4.6600E+07 

1.2000E-04 

3.1416E+00 

2.2982E+04 

4.4801E+07 

1.3000E-04 

3.3692E+00 

2.2542E+04 

4.3104E+07 

1.4000E-04 

3.5925E+00 

2.2119E+04 

4.1502E+07 

1.5000E-04 

3.8116E+00 

2.1712E+04 

3.9987E+07 

1 . 6000E-04 

4.0268E+00 

2.1319E+04 

3.8554E+07 

1.7000E-04 

4.2381E+00 

2.0940E+04 

3. 719GE+07 

1.8000E-04 

4.4456E+00 

2.0575E+04 

3.5909E+07 

1.9000E-04 

4.6496E+00 

2.0222E+04 

3.4688E+07 

jj  - 

2 . OOOOE-04 

4 . 8501E+00 

1.9881E+04 

3.3528E+07 

2.1000E-04 

5.0473E+00 

1.9551E+04 

3.2423E+07 

2.2000E-04 

5.2412E+00 

1.9232E+04 

3.137GE+07 

2.3000E-04 

5.4320E+00 

1.R924E+04 

3.0377E+07 

% 

■ 
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Table  A-V.  Kinematic  Results  (B>0, 


A=C=0) 


TEST  PRO]’, LEM  3 


MASS-  8. 

000E+00  RADIUS- 

6.000E-01  i:;i7 

VEL-  3.000E+04 

A-  0. 

000E+00  B- 

1.800E+03  C- 

0.000E+0O 

A3-  O.OOOE+OO  B3- 

KINEMATIC  RESULTS 

2.544E+03  C3- 

O.OOOE+OO 

TrtE(SEC) 

PE1I  (at) 

VEL(CM/SEC) 

del  E/DEL  s 

O.OOOOE+OO 

O.OOOOE+OO 

3.0000E+04 

7.4340E+07 

1.0000E-05 

2.9621E-01 

2.9248E+04 

7.4422EJ07 

2.0000E-05 

5.0498E-01 

2 . 8511E+04 

7.2532E-*-07 

3.0000E-05 

3.6650E-01 

2.7795E+04 

7.0729E+07 

4 . 0000E-O5 

1.1409E+0P 

2.709CE+04 

6.8952E+07 

5.0000E-O5 

1.4084E+00 

2.6415E+04 

6.7220E+07 

f..0000E-05 

1.6693E+00 

2.5752E+04 

6.5531E+07 

7.0000E-05 

1.9235E+00 

2.5105E+04 

6. 3884E+07 

3.0000E-05 

2.1714E+00 

2.4474E+04 

6.2279E+07 

9.0000E-05 

2.4131E+00 

2.3859E+04 

6.0714E+07 

1.0000E-04 

2.6487E+00 

2.3259E+04 

5.9109E+07 

1.1000E-04 

2.8783E+0D 

2.2675E+04 

5. 7701E+07 

1.2000E-04 

3.1022E+00 

2.2105E+04 

5.6252E+07 

1.3000E-04 

3. 3205E+00 

2.1550E+04 

5.4838E+07 

1.4000E-04 

3. 5333E+00 

2.1000E+04 

5 . 3460E+07 

1.5000E-04 

3.7407E+00 

2.0480E+04 

5.2117E+07 

1.6000E-04 

.3.9429E+00 

1.9966E+04 

5.0808E+07 

1.7000E-04  . 

4.1401E+00 

1.9464E+04 

4.9531E+07 

1.3000E-O4 

4.3323E+00 

1.8975E+04 

4.8286E+07 

1.9000E-04 

4.5196E+00 

1.8498E+04 

4.7073E+07 

2 . 0000E-04 

4.7023E+00 

1.8034E+04 

4.. 38908+07 

2.1000E-04 

4.8804E+00 

1.7580E+04 

4.4737E+07 

2.2000E-04 

5.0539E+00 

1.7139E+04 

4.3613E+07 

2.3000E-04 

5.2232E+00 

1.6708E+O4 

4.2518E+07 
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Table  A-VI.  Kinematic  Results  (C>0,  A=B=0) 


TEST  PROBLEM  4 


MASS**  3. 

0008+00  RAPITTS- 

G.900E-01  I! ’IT 

VKL=  3. 0008+04 

A=  0. 

OOOE+O''  R= 

O.OOOE+OO  C- 

5. 3998+08 

A>  0.0008+00  B3- 

KINEMATIC  RESULTS 

0.0008+00  C3= 

7. 6338+07 

ti:ie(sec) 

PEN  (Cl) 

vtl(cm/sec) 

PEL  E/7) El.  S 

o.oooon+on 

O.OOOOE+OO 

3.0000E+04 

7. 63378+07 

1.00008-05 

2.96138-01 

2. 92388+04 

7.60378*0? 

2.00008-05 

5. 047 38-01 

2.84738+04 

7. 6+ 378+07 

3.00008-05 

o.5564e-oi 

2.7709E+04 

7.63378+07 

4.0000E-05 

1.1339E+00 

2.8946E+04 

7.63378+07 

5.00008-05 

1.40458+00 

2.61338+04 

7. 63378+07 

■'.00008-05 

1.6625E+00 

2.54198+04 

7. 63378+07 

7.00008-03 

1.9129E+00 

2. 46568+04 

7.63378+07 

8.00008-05 

2.1557E+00 

2. 3892E+04 

7.63378*07 

9. OOOOE-05 

2.3908E+00 

2. 312PE+04 

7.63378+07 

1. 00008-^4 

2.61838+00 

2.2366E+04 

7.63378+07 

1.10008-04 

2.R3R18+00 

2.16028+04 

7.63378+07 

1.20008-04 

3. 05038+00 

2.O339E+04 

7. 6 3378+0? 

1.30908-04 

3.25498+O0 

2.0076E+04 

7.63378-H)7 

1.4000E-O4 

3.4518E+00 

1.93128+04 

7.63378+07 

1.50008-04 

3.64118+00 

1.85408+04 

7.6337E+07 

1.6900E-04 

3.82288+00 

1.77358+04 

7.6337r+07 

1.70OQE-O4 

3. 99698+00 

1.7022E+04 

7.63378*07 

1.R000E-04 

4.18338+00 

1. 62598+04 

7.63378+07 

1.90008-04 

4. 3221K+0Q 

1. 5495E+04 

7.63378*07 

7..0000E-04 

4.4732E+00 

1.4 7 328+04 

7.63378+07 

2.10005-04 

4.6167E+00 

1.39608+04 

7.63378+07 

2.2000E-O4 

4.75288+00 

1. 32058+04 

7.63378*07 

2.30008-04 

4.8808E+00 

1.24428+04 

7.63378*07 

j 
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PENETRATION  (CENTIMETRES) 


APPENDIX  B 


NEWT  PROGRAM  LISTING 

4 


\ 


i 


1REM  NEWT3 

2REM  F=-A*V*V-B*V-C 

3 INPUT  "TITLE", Z2$,Z3$,Z4$,Z5$ 

4SELECT  LIST  211 

5SELECT  PRINT  211 

GSELECT  R 

lODiri  T1  (50)  ,V1  (50) , SI  (50) 

15 INPUT  "A=",A 

20 INPUT  "B=",B 

25 INPUT  "C=",C 

30 INPUT  "RADIUS3" ,R 

35 INPUT  "MASS=",M 

40 INPUT  "INITIAL  VELOCITY3", VO 

45 INPUT  "MO.  TIMES3", N1 

50INPUT  "1  FOR  READ-IN  TIMES", F2 

55 IF  F2=0THEN  75 

60F0R  1=1 TO  N1 

65 INPUT  "NEXT  TIME3",T1(I) 

70UEXT  I 

75A3=*PI*0.5*R*R*A/M 
30E  3*6 . 0*rP I*R*B/M 
35C3=#PI*R*R*C/M 
90 IF  F?]OTIIEN  125 
95T=0.0 

100INPUT  "TIME  INTERVAL3", Z1 

105F0R  1=1  TO  N1 

HOT!  (I)=T 

1 15T=T+Z1 

120NEXT  I 

125 IF  A]OTHEN  134 

127 IF  B]OTUEN  130 

12007=001:  GOTO  220 

130IF  C]OTHEN  132 

13107=010:  GOTO  220 

132Q7=011 : GOTO  220 

134 IF  P,]0THEN  138 

135 IF  C]OTUEN  137 

136Q7=100:  GOTO  154 
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139Q7=H0:  GOTO  154 
140Q7=111 

154REM  CALCULATE  B*B-4*A*C  FLAG 

155F1=B3*B3-4.0*A3*C3 

160IF  FI [0. OTHER  175 

1C5IF  FI  JO. OTHER  130 

170F1=3:G0T0  220 

175F1=2:G0T0  220 

luOFl=l 

220REI1  PRINT  INPUT  AND  CAPTIOUS 
222PRIHT  Z2$,Z3$,7.4$,Z5$ 

223Z2$="MA5S=" 

224Z3$="RAUIUS=" 

225Z4$="IIIIT  VEL=" 

220PR I RTUS I f IG  227 , Z25 , M , Z3$ , R , Z4$ , VO 

r.r*##!!!!  mmm  ?.###!! n mmm  t.wriiii 

223Z2$="A=" 

229Z3$»"B»" 

230745=' "C=" 

231  PRI'ITUSiriG  227,Z25,A,Z3$,B,7.4$,C 
232Z25="A3=" 

333Z3$="B3=" 

734745= "C 3= " 

235PRIRTUSIRG  227,Z2$,A3,Z3$,B3,Z4$,C3 
245GOSUR  Ml 

246 PRINT  "KINEftATIC  RESULTS" 

250Z1 5=  "TIME  ( SEC) " : Z2$="PEf!(Cfl) " : Z3$="VEL  (Cil/SEC ) " 

251Z4$=  "UCL  F./9EL  S" 


’’55PRI  f ITUS  I MG  260 , Z1  $ .Z2$ , 7.3$  ,Z4$ 


2w%mmm  mm* 

265GGSIIB  Ml 
270F0R  1=1  TO  111 
275IF  Q7]1THEII  795 
280G0SIJB  'IKTI(I)) 
2C5G0SIJC  '21 
290G0T0  340 
295 IF  Q7]11TIIEIJ  315 
3O0G0SUB  * 1 3 (Tl (I)) 
3O5G0SIJB  '23 
31OG0TO  340 
31  GIF  Q7]lllTHEtl  335 
320G0SIJB  *ir,(Tl  (I) ) 
325G0SIIB  '26 
330G0T0  340 
335A9=2:G0T0  430 
340HF.XT  I 

342 IF  Tl (1 )]0.0THEN  345 
343S1 (1 )=0.0 
345F0R  1=1  TO  111 
346 IF  VI  (I)[0TI!ER  352 
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350PRINTIISING  355.T1 ( I) ,S1 ( I) ,V1 ( I ) ,A3*V1 ( I ) !2+B3*Vl ( I )+C3 
351  GOTO  3f»0 

352PRI rfT  "NEGATIVE  VELOCITY" 

353N1 *1—1 
354G0T0  365 

355*#.####! II!  #•####!!!!'  #.####111!  #.####! ! I ! 

360NEXT  I 

365G0SUB  '01 

367 FOR  1=1 TO  5 

368PRINT 

369NEXT  I 

370G0SUB  '01 

371Z2$S"A3*V*V":  Z3$="B3*V":  Z4$="C3" 

372PR1NT  "DIFFERENTIAL  ENERGY  LOSS (ERGS/CM)" 

373PRI NTUS ING  374 , Z1 $ , Z2$ , Z3$ , Z4$ 

374  *#########  mm  ####  ## 

375G0SUB  '01 
376F0R  1=1 TO  N1 

378PRINTUSING  355 ,T1(I) .A3*V1 ( I )*V1 ( I ) ,B3*V1 (I),C3 
380NEXT  I 
382G0SUB  '01 


334F0R  1=1  TO  5 
336PRINT 
383 NEXT  I 
389G0SUB  '01 

390 PR I NT  "CUMULATIVE  ENERGY  LOSS (ERGS)" 
391Z2$="INT(A3*V*V)DS" 

3927.3$="INT(B3*V)DS" 

393Z4$="INT(C3)DS" 

396 PRINTUSING  393,Z1$,Z2$,Z3$,Z4$ 

398 %mmm  mtmmw  mmmn  mmm 

399G0SUB  '01 
400Z6=T1(1) 

402Z7=0.5*A3*(V0!2+V1 (1 ) !2)*S1 (1 )*M 

40478=0. 5*B3*(V0+V1 (1 ) )*S1  (1  )*M 

406Z9=C3*S1 (1 )*M 

408PR I NTUS I NG  355 , Z6 , Z7 , Z3 , Z9 

410F0R  I=2T0  N1 

412ZG«T1 (I) 

414Z7=Z7-K).5*A3*(V1 ( 1-1 ) S2+V1 ( I ) 12)*(S1 ( I )-Sl ( 1-1 ) )*M 
41 6Z8=Z3+0.5*B3*(V1 ( 1-1 )+Vl (I) )*(S1 (I)-Sl ( 1-1 ) )*M 
413Z9=Z9+C3*(S1 (I)-Sl ( 1-1 ) )*M 
420PRINTIISING  355,Z6,Z7fZ8,Z9 


f. 


422NEXT  I 
424GOSUD  '01 


430END 

51 000EFFJI  'll  (T) 

5102REII  CALCULATES  VELOCITY  07=001 

5105V=V0-C3*T 

5110V1 ( I)=V 


i 


j 

i 


i 

J 


I 

? 

I 
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5115RETURN 

5300DEFFN'13(T) 

5302REM  CALCULATES  VELOCITY  Q7-011 
5305V»B3*T  :V*EXP(-V) 
5310V=V*(B3*V0+C3) 

531 5V=(V-C3)/B3 
5320V1 (I)=V 
5325RETURN 
56U0DEFFN'1G(T) 

5C02REM  CALCULATES  VELOCITY  Q7=lll 

5G05IF  FI ]1 THEN  5645 

5G07REM  LOO  FORM 

561  OR! »D3! 2-4 .0*A3*C3 : R1 =SQR ( R1 ) 

5G15F3=2.0*A3*V0+R3+R1 

5G20F3=F3/(F3-2.0*R1) 

5625V=EXP(R1*T):  V=V*F3*(R1-B3) 

5630V=V+R3+R1 

5G35  V=V/(2.0*/'.3*(F3*EXP(Rl*T)-l  .0) ) 

5640G0T0  5G94 

GG45IF  FIFTEEN  5GC0 

5C47RHI  TRIG  FORM 

j65( )R  1 =/■■ . 0*A3*C3-R 3 ! 2 : R1-SQR(R1 ) 

5655 V=2 . H*A 3*V9+l  3 : V= V/R 1 

5660V® A!  CTAfi(V) : V-V-Rl*T/2.n 

!>6G5V=TAh(V)  -B3/R1 

5670Vb0. 5*R1 *V/A3 

5G75G0T0  5G94 

5600 IF  FI  ]3TIIEN  5G92 

5632RC1  ALGEBRAIC  FORM 

5605  V -2*A 3*V  0+5  3 

66S7V»2*A3*( 2+T*V ) 

5600V-  ( 1-*A3*V.>R3*T*  ( r.*A3*V0+H3 ) ) /V 

5692A'>l:GGTi.  5696 

3694V 1 ( I )=V 

GCOGRETilRN 

6100  nCFFN'21 

6132i:C;  CALCULATES  PENETRATIO!!  07=001 

610GS=V0*V0-V*V 

6110S=0.5*S/C3 

G115S1(I)=S 

6120RETURN 

G33..iirrFN'23 

6302RE.  i CALCULATES  PENETRATION  Q7-011 

6305S«P.3*V0+C3 

631  viS=  ( C 3*V+C3  )/S:S«LOG(S) 

G3155=C3*S/B3 

G320S=(V0-V+S)/B3 

6325S1 (I)=5 

G330RETURN 

6600DEFFfl'26 
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6602REM  CALCULATES  PENETRATION  07-111 
6625Z1 -A3*V*V+D3*V+C3 
6630Z1 » ( A3*V0*V0+B3*V0+C3)/Z1 
6635Z1«L0G(Z1) 

6637S«B3*T1(I) 

6640S«0.5*(Z1-S)/A3 
6740  S1(I)-S 
C745RETURN 
7000DEFFN'01 

7002REM  PRINTS  DASHED  LINE 

70057.9$-" " 

7010PRINT  Z9$,Z9$,Z9$,Z9$ 

701 5RETURN 
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